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Abstract. Using the IRAM interferometer we have detected J=4^3 and 5^4 CO line emission toward the radio 
quiet quasar PSS 2322+1944. At a redshift of zco = 4.1199 this is the fourth and strongest detection of CO at 
z > 4. The velocity-integrated CO J=4^3 and J=5^4 line fluxes are 4.21 ± 0.40 and 3.74 ± 0.56 Jy kms -1 , and 
the linewidth is « 300 kms -1 . The CO J=10^9 was searched for but not detected with an upper intensity limit 
of 30 mJy. The 1.35 mm (250 /im rest wavelength) continuum flux density is 7.5+1.3 mJy, in agreement with 
previous bolometer measurements at 1.2 mm with the 30-m IRAM telescope. The 3 mm (580 /im rest wavelength) 
continuum is not detected with a 3 a upper limit of 0.7 mJy. We also report observations of the 450 /im continuum 
in PSS 2322+1944 using the SCUBA array at the JCMT. The quasar was detected with a 450 /im flux density 
of 79 ± 19mJy. At the angular resolution of 4" 8 x 2"1 at 1.3 mm and 6" 2 x 4" 9 at 3.2 mm, the interferometer 
observations do not show evidence of spatial extension in the continuum or CO line emission. Assuming no 
gravitational magnification, we estimate a molecular gas mass of « 2.5 x 10 11 M . The molecular gas is warm 
(40 - 100 K) and dense (10 3 5 - 10 41 cm" 3 ). The infrared- to- CO luminosity ratio is « 185 L (Kkms -1 pc 2 ) -1 , 
comparable to the values found for ultraluminous infrared galaxies. The detection of CO emission in this high 
redshift QSO provides further evidence that the radio emission and the millimeter to submillimeter continuum 
emission are predominantly powered by a starburst which is coeval with the AGN activity. 

Key words. Galaxies: formation - Quasars: emission lines - Quasars: individual: PSS 2322+1944 - Cosmology: 
observations - Cosmology: early Universe - Radio lines: galaxies 



1. Introduction 

The study of dust and molecular gas in sources at high 
redshift has opened up new ways to probe the physical 
conditions during the early evolution of galaxies and to 
study the star formation history in the early universe. 
Spatially resolved observations of the molecular gas con- 
tent in the host galaxies of high-z quasars are a key 
for understanding the relationship between black hole 
formation and spheroidal galaxy formation. At redshifts 
0.04 < z < 0.27, the CO J=1^0 line has been detected 
in 14 QSOs (Evans et al. 2001; Casoli & Loinard 2002) 
indicating large masses of molecular gas (a few 1O 9 M ) 
fueling both the AGN and star formation in the QSO host 
galaxies. At higher redshifts, although the lines are broad 
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and weak and precise redshifts are difficult to predict, in 
the last decade 12 CO emission was detected in 15 galax- 
ies at 1.44 < z < 4.69 (see Guilloteau 2002 and reference 
therein; Andreani et al. 2000; Barvainis et al. 2002). Most 
of the detections were obtained for high rotational transi- 
tions of CO (J=3^2 to 9^8) lines which are redshifted 
to millimeter wavelengths. In a few cases, the J=2^1 and 
1— >0 lines were detected at centimeter wavelengths (Carilli 
et al. 1999, 2002; Papadopoulos et al. 2001). These CO 
observations imply that the high redshift galaxies contain 
large amounts, a few 10 10 M , of molecular gas which are 
predominantly excited by giant starbursts with star for- 
mation rates ~ lO 3 M0yr -1 (see, e.g., Guilloteau et al. 
1997, 1999). 

In recent years we have concentrated our search for 
high-z CO emission to QSOs with strong (sub) millimeter 
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dust continuum emission. Based on the first surveys of 
(sub) millimeter emission of high-z QSOs (McMahon et 
al. 1994; Omont et al. 1996a), we reported three detec- 
tions of CO emission in z > 4 sources: BR 1202—0725 
(Omont et al. 1996b), BRI 1335-0415 (Guilloteau et al. 
1997) and BRI 0952-0115 (Guilloteau et al. 1999). We 
recently extended these surveys for continuum emission 
using the Max-Planck millimeter bolometer (MAMBO) 
array at the 30-meter (Omont et al. 2001; Carilli et al. 
2001a) and SCUBA at the JCMT (Isaak et al. 2002). 
Continuum emission was detected toward 35 new high- 
z QSOs: Omont et al. (2001) and Isaak et al. (2002) 
observed a total of 75 bright quasars at z > 3.8 from 
the Palomar Sky Survey (PSS) sample and detected 19 
sources; Carilli et al. (2001a) made a deeper survey of a 
sample of 41 QSOs at z > 3.7 selected from the Sloan 
Digital Sky Survey (SDSS) and detected 16 of them. 

PSS 2322+1944, discovered by Djorgovski et al. (in 
prep.), is the brightest QSO in the survey of z > 3.8 PSS 
quasars by Omont et al. (2001) and Isaak et al. (2002), 
with flux densities of 9.6 ± 0.5 and 22.5 ± 2.5 mJy at 
250 and 350 GHz, respectively. The luminosity and im- 
plied dust mass are estimated]] to be ~ 2.7 x 10 13 L 
and ~ 1.6 x 10 9 M (see Sect. 3). Under the assump- 
tion that the far-infrared luminosity arises exclusively 
from young stars, the implied star formation rate is a few 
10 3 M o yr- 1 . PSS 2322+1944 was also detected with the 
VLA at 1.4 GHz with a total flux density of 98 ± 20/iJy 
being spatially extended on a scale of 1.5" (Carilli et 
al. 2001b). The ratio of radio to millimeter flux agrees 
with the radio-to-far infrared correlation for star-forming 
galaxies. PSS 2322+1944 is also an exceptional optical 
source as described by Djorgovski et al. (in prep.). From 
a Keck optical spectrum a redshift was derived from mul- 
tiple, associated narrow absorption lines, yielding z opt = 
4.11075 ±0.0005. 



2. Observations and Results 

Observations were made with the IRAM Plateau de Bure 
interferometer during 5 nights in June 2001. We used the 
standard CD configuration (4 antennas) which results in a 
beam of 4'/ 8 x 2'/l at 1.3 mm and 6'/ 2 x 4'/9 at 3.2 mm. Dual 
frequency receivers were used to search simultaneously for 
CO J=4^3 and dust emission at 3 and 1.35 mm, and in 
a second series of observations for CO J=5^4 and 10^9 
emission. During the first observations, the receivers were 
tuned at 90.03 and 90.06 GHz to search for the J=4^3 
emission line. A subsequent series of observations was done 
with the 3 mm receiver centered at 112.6 GHz (at the 
position of the redshifted CO J=5^4) and the 1.3 mm 
receiver tuned at 225.0 GHz to search for CO J=10^9. 
The 3 mm receivers were tuned in single sideband and the 
1.3 mm receivers in double sideband. Typical SSB sys- 



tem temperatures were ~ 150 K at 3 mm and ~ 400 K 
at 1.3 mm. Amplitude and phase calibration were done 
using 3C 454.3 (6.1, 5.3 Jy) and 2230+114 (2.1, 1.2 Jy) 
where the number in round brackets are the flux densities 
at 3 and 1.3 mm, respectively, at the date of the obser- 
vations. Phase noise was stable on all baselines both at 
3 mm (rms 5°-17°) and at 1.3 mm (rms 13°-45°). During 
the observations, the water vapour varied between 3 and 
6 mm and the seeing conditions varied between 0'/8 and 
1V3, i.e. much smaller than the synthesized beams at 3 mm 
and 1 mm. The total integration time was about 18 hours 
for the CO J=4^3 line, and 12 hours for the CO J=5^4, 
the CO=10^9 and the 1.35 mm continuum. The data 
were reduced, calibrated and analyzed using the standard 
IRAM programs CLIC and MAPPING. The final spectra 
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Fig. 1. Spectra of the CO J=4^3 and J=5^4 lines to- 
ward PSS 2322+1944. Error bars are ±la. The velocity 
scale corresponds to frequencies of 90.05 and 112.55 GHz, 
or a redshift z C o = 4.1199 ± 0.0008. The redshift de- 
rived from multiple, associated narrow absorption lines 
in a Keck optical spectrum is z op t = 4.11075, which is 
indicated by an arrow. The weak continuum was not sub- 
tracted. 
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Fig. 2. (a) Channel maps of the CO J=4^3 line toward PSS 2322+1944. The channel width is 35 MHz, i.e., 120 km s -1 . 
The contour step is 1.33 mJy beam -1 . The offset velocity is indicated in the upper left corner of each map. The zero 
velocity corresponds to 90.05 GHz, i.e., z = 4.1199. The beam of 6.25 x 4.98 arcsec 2 is shown in the lower left panel. 
(b) The lower right panel shows the 1.35 mm continuum map. Contour steps are 2 mJy beam -1 , and the r.m.s. noise 
1.5 mJy beam -1 . The beam of 4.85 x 2.14 arcsec 2 is shown in the lower left corner. Offsets are relative to the optical 
(and radio) position (indicated by a cross), R.A. 23:22:07.25, Dec. 19:44:22.08 (J2000.0). 



are shown in Fig. 1 and, for the CO J=4^3 emission, 
Fig. 2a also shows the velocity channel maps. 

PSS 2322+1944 was detected in the continuum at 
1.35 mm with a flux density of 7.5±1.3 mJy, as derived 
from the image sideband measurements. This value is con- 
sistent with the flux density measured at 1.20 mm with 
MAMBO at the 30-meter (9.6±0.5 mJy) when the spectral 
index of dust emission is taken into account. Within the 
astrometric uncertainties (±0'/3), the continuum source 
corresponds is to the optical (and radio) position, RA 
23:22:07.25, Dec 19:44:22.08 (J2000.0) from Carilli et al. 
(2001b). At the angular resolution of our observations, 
the source is not resolved (Fig. 2b). PSS 2322+1944 is not 
detected in the continuum at 3 mm with a flux density 
of 0.40 ± 0.25 mJy as derived from the line-free channels. 
Using a dust spectral index of 3.5, this is consistent with 
the 3 mm continuum value of « 0.3 mJy expected from the 
1.35 mm flux density (Fig. 3). 

PSS 2322+1944 was also observed with the wide-band 
450 /im and 850 /im filters on SCUBA at the JCMT 
in December 2001 under good and stable weather con- 
ditions with T850 /2m ~ 0.06. The observations, calibration 
and data reduction were done as explained in Isaak et 
al. (2002). The source was detected in the continuum at 
450 jam (corresponding to a rest wavelength of 88 /im) 
with a flux density of 75±19 mJy. At 850 /im, the flux 
density of PSS 2322+1944 is 24±2 mJy, which is entirely 
consistent with the 22.5±2.5 mJy value reported in Isaak 



et al. (2002). Fig. 3 presents the spectral energy distri- 
bution of PSS 2322+1944 with all the currently available 
photometric data. 

Both the CO J=4^3 (^ rest = 461.0408 GHz) and 
J=5^4 (^ r est = 576.2679 GHz) are clearly detected to- 
ward the position of the continuum emission. The J=4^3 
and 5^4 lines are found at frequencies of 90.05 and 
112.55 GHz (Table 1) corresponding to a redshift of 
zco = 4.1199 ± 0.0008, close to the redshift derived from 
the optical spectrum (z opt = 4.11075). The difference cor- 
responds to a velocity difference of 530 km s -1 (Fig. 1). 
The CO J=10— >9 emission line was not detected with a 
3<j upper intensity limit of 30 mJy (Table 1). The inte- 
grated line fluxes of the CO J=4^3 and J=5^4 lines 
are 4.24 ± 0.33Jykms -1 and 3.74 ± 0.56, respectively. 
PSS 2322+1944 is the strongest CO emitter at high z, 
even stronger than the lensed quasar APM 08279+5255 
(Downes et al. 1999) - see Table 2. From Gaussian fits, 
the line widths are found to be « 300 km s -1 (Table 1) 
which is comparable to the width in the three z > 4 CO 
sources detected so far (Table 2). 

The channel maps in the CO J=4^3 emission line 
(Fig. 2a) (as well as the J=5^4 channel maps not shown) 
do not show evidence for extension and/or position shifts 
with velocity. The 3.2 mm CO data are limited by the 
angular resolution of the observations (6'/ 2 x 4'/ 9) and 
other array configurations will be needed to probe further 
whether the emission of PSS 2322+1944 is extended. 
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Table 1. Observed Properties of CO Lines in PSS 2322+1944 



Line 


^obs 


Peak Int. 


AVfwhm 


^co 


Continuum 








[GHz] 


[mjy] 


[kms- 1 ] 


[Jy km s _1 ] 


[mjy] 


[10 11 Kkms- 1 pc 2 ] 


[1O 8 L ] 


CO(4^3) 


90.05 


10.5 


375+41 


4.21+0.40 


0.40+0.25 


2.0 


6.3 


CO(5^4) 


112.55 


12.0 


273+50 


3.74+0.56 




1.1 


7.0 


CO(10^9) 


225.00 


< 30 




<5.2f 


7.5 + 1.3 


< 0.4 


< 19 



NOTE. - 1 Adopting a line width of 300 kms -1 

Table 2. Comparison of the CO Results of PSS 2322+1944 with other z > 3.5 Sources 



Source 


z 


CO(4 


-3) 


CO(5 


-4) 


Ref. 






^co 


AVfwhm 


^co 


AVfwhm 








[Jy km s _1 ] 


[km s" 1 ] 


[Jy km s _1 ] 


[km s" 1 ] 




PSS2322+1944 


4.12 


4.21+0.40 


375+41 


3.74+0.56 


273±50 


[1] 


BRI 1335-0415 


4.41 






2.80+0.30 


420±60 


[2] 


BRI 0952-0115 


4.43 






0.91+0.11 


230±30 


[3] 


BR 1202-0725 


4.69 


1.50+0.30 


280+30 


2.40+0.30 


320±35 


[4,5] 


APM 08279+5255 


3.87 


3.70+0.50 


480+35 






[6] 


4C 60.07 


3.79 


2.50+0.43 


> 1000 






[7] 


6C 1909+722 


3.53 


1.62+0.30 


530+70 






[7] 



NOTE. - [1] This paper [2] Guilloteau et al. (1997) [3] Guilloteau et al. (1999) [4] Ohta et al. 1996 [5] Omont et al. (1996b) [6] 
Downes et al. (1999) [7] Papadopoulos et al. (2000) 



3. Discussion 

The observed CO J=4^3 and 5^4 line fluxes imply 
intrinsic CO line luminosities, L' co = 2.0 and 1.1 x 
10 11 K kms -1 pc 2 , respectively, or expressed in solar lu- 
minosities, L' co = 6.3 and 7.0 x 10 8 L (Table 1). Lqq is 
proportional to line brightness (Rayleigh- Jeans) temper- 
ature integrated over the area of the source: the ratio of 
luminosities in two CO transitions originating from the 
same area is a measure of the line brightness ratio and 
therefore an indicator of the physical conditions in the 
molecular gas (see Solomon et al. [1997] for a discussion 
of L' co and L co ). 

From the L' co we derive brightness (Rayleigh- Jeans) 
temperature ratios T 6 [CO(4^3)]/T 6 [CO(5^4)] = 1.76 ± 
0.25 and T 6 [CO(10^9)]/T 6 [CO(5^4)]<0.34. The former 
ratio is comparable to the values of 1.38 derived for the 
z = 2.56 quasar H 1413 (the Cloverleaf) - see Barvainis 
et al. (1997) - and of 1.02 found for the z = 4.69 quasar 
BR 1202-0725 (Omont et al. 1996a). These ratios indicate 
that the molecular gas in PSS 2322+1944 must be warm 
and dense. Using the results of a Large Velocity Gradient 
model for high z galaxies (Sakamoto 1999), we used the 
brightness temperature ratios to constrain the molecular 
gas density in PSS 2322+1944. For a gas kinetic tempera- 
ture of « 40 — 100 K, in accordance with the temperature 
derived from the thermal dust emission spectrum (see be- 
low), and adopting a value of 10 -6 pc/(kms _1 ) for the 
CO abundance per unit velocity gradient of the molecu- 
lar gas (X (CO) /dv/dr), the T 6 [CO(4^3)]/T 6 [CO(5^4)] 
ratio constrains the density of the molecular gas to be of 
the order of 10 3,5 — 10 4 1 cm -3 . These values for the gas 



temperature and density are also consistent with the line 
brightness ratio T 6 [CO(4^3)]/T 6 [CO(1^0)] « 1.4 derived 
from the CO (1— »0) line flux measured with the VLA in 
PSS 2322+1044 (Carilli et al. [in prep.] - a more detailed 
analysis of the excitation conditions will be given in this 
paper). The gas density in PSS 2322+1944 is compara- 
ble to the high densities (n(H2) ~ 10 4 cm -3 ) found in the 
highly excited gas in the nuclear region of the starburst 
galaxy M 82 (Glisten et al. 1993; Mao et al. 2000) and in 
nearby ultraluminous infrared galaxies (e.g., Solomon et 
al. [1992]). 

Assuming a standard Galactic value of 
4.6 M (K kms -1 pc 2 ) -1 for the conversion factor of 
molecular mass to -^co(i-o)' anc ^ usm & the above line 
brightness ratio T 6 [CO(4^3)]/T 6 [CO(1^0)], we infer a 
molecular gas mass of « 6 x 10 11 M . Based on a study 
of ultraluminous infrared galaxies (ULIRGs), Solomon 
et al. (1997) showed that the molecular gas mass to CO 
luminosity conversion factor is likely to be lower (by a 
factor & 3) in the extreme environments of luminous 
infrared galaxies where the molecular clouds are expected 
to be hotter and denser than in the Galaxy (see also 
Combes et al. 1999). In the case of Arp 220, Scoville et 
al. (1997) derived a conversion factor « 0.45 times the 
Galactic value. The above gas mass is therefore likely to 
be overestimated by a factor 2 to 3. In the following we 
will adopt a conversion factor of 1.8 M (Kkms -1 pc 2 ) -1 
(i.e., a factor of 2.5 times smaller than the Galactic 
conversion factor), which translates into a molecular gas 
mass of M H2 « 2.5 x 10 n M o for PSS 2322+1944. 

Compared with the dust mass Md us t ~ 1.6 x 10 9 M 
derived from the 1.35 mm continuum flux density (see 
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Fig. 3. Spectral energy distribution of PSS 2322+1944. 
Triangles: 850 /im (Isaak et al. 2001); 1.2 mm (Omont 
et al. 2001); 1.4 and upper limit at 5 GHz (Carilli et 
al. 2001b); squares: 450 /im, 1.35 mm and the 3a upper 
limit at 3 mm (this paper). For comparison, the radio- 
to-infrared spectral energy distribution of the starburst 
galaxy M 82 is shown, red-shifted to z = 4.12 and normal- 
ized to the flux density of PSS 2322+1944 at the observed 
wavelength of 850 jam: the crosses show all the currently 
available photometric data, and the continuous line repre- 
sents the ISO LWS spectrum (from Colbert et al. 1999). 
The left- and right-hand flux density scale are adapted for 
PSS 2322+1944 and M 82, respectively. 

Omont et al. 2001), the implied gas to dust mass ratio 
for PSS 2322+1944 is M(H 2 )/M dust « 150, comparable 
to the ratios derived in bright nearby galaxies (Dunne et 
al. 2000), ULIRGs (Solomon et al. 1997), or other high 
redshift CO sources (Guilloteau et al. 1999). This similar- 
ity is also reflected in the Lfir/£q O (i-0) ra ti° s °f these 
sources. In the case of PSS 2322+1944, Lfir/^co(i-o) ~ 
185 L© (K kms -1 pc 2 ) -1 , close to the median value of 160 
found for a sample of 37 ULIRGs by Solomon et al. (1997). 

Figure || presents the radio to infrared flux densities 
of PSS 2322+1944 (including the new measurement at 
450 /im and 1.35 mm) compared with the spectral energy 
distribution of M 82, which we red-shifted to z = 4.12 
and normalized to the flux density of PSS 2322+1944 at 



the observed wavelength of 850 /im. Note that M 82 is 
about 100 times less luminous than PSS 2322+1944. The 
two spectral energy distributions are remarkably similar 
in the far-infrared/radio, although the peak intensity of 
PSS 2322+1944 is somewhat higher. For a dust emissiv- 
ity index (3 = 1.6, the dust emission is well fitted with 
a temperature of « 47 K, which is similar to the dust 
temperature derived in M 82 under the same assumptions, 
i.e. Td w 45 K. This further supports the idea that in 
PSS 2322+1944 the starburst dominates the dust heat- 
ing and the radio emission (see also Omont et al. [2001] 
and Carilli et al. [2001b]). Further confirmation may come 
from a comparison of their mid- and near-infrared spectra. 
Similarities with the far-infrared and radio SED of M 82 
were also noted for the QSOs BRI 1335-0417 at z = 4.4 
(Carilli et al. 1999) and BR 1202-0725 at z = 4.7 (Kawabe 
et al. 1999; Yun et al. 2000). 

Djorgovski et al. (in prep.) found that the quasar 
PSS 2322+1944 has a close companion (separated by 
about 1") with the same spectrum within the measure- 
ment uncertainties. It is not yet known if this is a case 
of a gravitational lensing, or a binary quasar. In the for- 
mer case, the luminosities and dust masses deduced here 
should be lowered by a factor of a few. In the latter case, 
we may be seeing an interaction or a merger of two quasar 
hosts. 

4. Conclusions 

With an apparent CO luminosity greater than that of the 
strongly lensed quasar APM 08279+5255, the z = 4.12 
quasar PSS 2322+1944 is the strongest CO emitter de- 
tected to date at redshifts larger than 3.5. Assuming no 
gravitational magnification, we estimate a molecular gas 
mass of « 2.5 x 10 11 M©, and a far-infrared luminosity of 
« 2.7 x 10 13 L©. The spectral energy distribution and large 
luminosity suggest that a massive starburst takes place in 
PSS 2322+1944, which may be related to the formation 
of the core of an elliptical galaxy. 

The exceptional brightness of PSS 2322+1944 makes 
it a good target for further observations of other CO tran- 
sitions, in particular those from lower levels which will 
constrain the physical conditions of the bulk of the molec- 
ular gas. Higher spatial resolution measurements are also 
needed to see whether the line or continuum emission is 
extended, as it was seen for the radio and optical emission 
(Carilli et al. 2001b; Djorgovski et al., in prep.). 

The detection of CO in another strong (sub) millimeter 
continuum high-z quasar confirms a frequent correlation 
between the 3 mm CO peak intensity and the 1.3 mm 
continuum flux, with a typical ratio of about unity. This 
relation shows that systematic searches for CO emission 
in strong thermal dust continuum quasars are promis- 
ing with current instrumentation, provided that the red- 
shift is known with high enough accuracy. Systematic 
(sub) millimeter continuum surveys of high redshift, radio- 
quiet quasars are therefore needed to find strong contin- 
uum sources towards which CO emission can be searched. 
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Such studies promise to further our understanding of the 
physical and chemical properties of the most energetic 
sources in the early Universe. Inter ferometric observa- 
tions, especially with ALMA and EVLA, will eventually 
be able to show the spatial distribution of the molecular 
gas and its relation to the stars and ionized gas. 
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